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A movable inclusion in an elastic material oscillates as a rigid body with six degrees of freedom. Dis- 
placement/rotation and force/moment tensors which express the motion of the inclusion in terms of 
the displacement and force at arbitrary exterior points are introduced. Using reciprocity arguments 
two general identities are derived relating these tensors. Applications of the identities to spherical 
particles provide several new results, including simple expressions for the force and moment on the 
particle due to plane wave excitation. 
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I. INTRODUCTION 

Faxen relations are named after Hilding Faxen who 
derived several identities for calculating hydrodynamic 
forces and torques on particles in low Reynolds number 
flows, e.g. Ref. 1. As an example of a Faxen relation, or 
law, the force and torque on a rigid sphere of radius a 
moving with velocity Vq and spinning with angular ve- 
locity u) in an unbounded fluid of viscosity fi and velocity 
field v(r) in the absence of the sphere are [Eqs. (3-2.46) 
and (3-2.47) of Ref. 2] 

F = 67r/ia[v(0) - vq] + 7r^a^V^v(0), (la) 

T = 47r/xa3[curlv(0) - 2u;]. (lb) 

The reader will note that the identities have as a spe- 
cial case the classical Stokes drag law, but they include 
additional effects caused by spatially variable flow fields. 
These and other Faxen relations for non-spherical par- 
ticles are based upon general integral identities relating 
the force and torque on the particle to the external flow 
field^"* . Although Faxen relations are commonly used in 
hydrodynamics and microfluidics, they seem to be essen- 
tially unknown outside that subject area. For instance, I 
am aware of only one mention^ of a Faxen type relation in 
elasticity and that one was is in regards to elastostatics. 

The objective of this paper is to develop similar ideas in 
the context of elastodynamics and in the process demon- 
strate their utility and wide application. Using dynamic 
reciprocity, a set of relations are first derived between 
the velocity or force of a particle in a solid matrix and 
the displacement or force at a distant point in the solid. 
These equations include but go far beyond the notion of 
particle impedance, which relates the force on a parti- 
cle to its velocity. Numerous applications of the general 
relations are obtained by considering spherical particles. 
Faxen-like relations are derived for the force and mo- 
ment on a spherical particle caused by plane wave inci- 
dence. Like their hydrodynamic counterparts, the elas- 
todynamic Faxen relations are simple in form. 
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The analysis here is the second in a series of papers de- 
veloping a simplifled algebra for calculating the radiation 
and scattering from inclusions in elastic and viscoclastic 
materials. In the previous paper^ the forced motion of 
a spherical particle in an clastic matrix was considered. 
Although this is a classical problem, originally solved by 
Oestreichcr^, it turns out that the dynamic impedance 
of the inclusion can be represented in a simplified man- 
ner. This was achieved^ through several lumped mass 
impedances for a spherical inclusion, in terms of which 
impedance or its inverse, admittance, has a simple form. 
The purpose of the present paper is to develop these ideas 
further, focusing on the interaction between the inclusion 
and remote points in the matrix. 



The plan of the paper is as follows. The dynamic prop- 
erties of an inclusion are defined in Section II, as well as 
some important quantities that are used throughout the 
paper: the displacement/rotation tensors U and W, the 
force/moment tensors €> and ^, and the impedances Z 
and Z. In Section III we prove the symmetry of these 
impedances, and derive two fundamental relations be- 
tween the displacement /rotation and force/moment ten- 
sors. The remainder of the paper focuses on the special 
case of spherical inclusions. The fundamental quantities 
for the spherical particle are presented in Section IV in a 
concise format using lumped mass impedances. Section 
V is the longest in the paper, as it contains numerous 
applications, discussion of limiting cases, and the new 
elastodynamic Faxen relations that are analogous to the 
classic hydrodynamic identities. The many results and 
their import are summarized in Section VI. 



Regarding notation, the time harmonic factor e~*"* 
is omitted but understood. Boldface quantities are ei- 
ther vectors or second order tensors. Vectors are usu- 
ally denoted by lower case, and tensors are capitalized, 
with the exceptions F and M which indicate force and 
moment vectors, respectively. The axial tensor axt(a) 
of the vector a is a skew symmetric tensor defined by 
axt(a)b = a A b. 
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II. INCLUSIONS AND RIGID BODY MOTION 

An inclusion in a solid matrix is defined to be the sur- 
face dVp of a finite volume Vp within which there could 
be a particle, or there could be some complicated "black 
box" with its own internal dynamics. The key feature 
of the inclusion is that its boundary dVp undergoes rigid 
body motion. In this sense the boundary is a rigid inter- 
face between the particle, whatever that may be, and the 
solid matrix. The term inclusion rather than particle is 
used throughout in order to remind us of this distinction. 

A. Tensor functions U, W, $ and "S* 

Rigid body motion has six degrees of freedom, which 
we characterize by two vector quantities: Up and Bp. Up 
is the rigid body displacement of the inclusion center of 
mass. Bp describes the rotation of the inclusion about the 
center of mass. The most general displacement possible 
for the inclusion is 

up = Up + 6>pAr, VredVp, (2) 

where r is the position relative to the center of mass. 
We will use the vector up to denote the total rigid body 
displacement, with Up reserved for the linear part (the 
entire development in this paper applies only to linear 
as distinct from nonlinear motion, so that the term lin- 
ear is synonymous with rectilinear). Note that Up has 
dimensions of length while Bp is dimensionless. For the 
sake of simplicity it is useful to consider the linear and 
rotational motions separately. Figure 1 shows the inclu- 
sion oscillating back and forth with linear displacement 
Up. In the absence of other sources of vibrational energy, 
the inclusion motion induces motion at every point r in 
the exterior region V = M.^ /Vp according to u = U(r)up, 
as depicted in Fig. 1. Here U is a second order tensor 
defined everywhere in the matrix. In the same way the 
particle displacement at r £ ^ caused by a pure rotation 
of the inclusion may be defined by a second order tensor 
W. In short, the tensors U and W relate the rigid body 
displacement of the inclusion to the displacement u(r) in 
the exterior solid medium V according to 

u(r) = U(r)up + W(r)0p. (3) 

We next define two dual tensor functions associated with 
force and moment, respectively. 

Consider the situation in which a point force of mag- 
nitude times direction equal to F acts a,t s € V, inducing 
motion of the inclusion. Figure 2. The tensors $ and ^ 
define the net force Fp and couple Mp on the inclusion 
caused by the point force according to 

Fp = *(s)F, Mp = *(s)F. (4) 



1. External impedances 

We introduce two impedance tensors: Z and Z, called 
external impedances because they depend upon the ex- 
terior properties of the solid matrix. 




U(r)Up 

FIG. 1. The inclusion undergoes time liarmonic linear dis- 
placement Up, resulting in displacement u = U(r)up at posi- 
tion r. 




FIG. 2. The time harmonic point force F is applied at s, 
resulting in the force Fp — $(s)F on the inclusion. 



The impedance Z relates the force acting on the in- 
clusion to the inclusion linear velocity, see Fig. 3 which 
is similar to the situation in Fig. 1. It is assumed that 
either force or velocity is controlled and the other is the 
dependent variable, and that there is no other excitation 
from sources in V. Thus, let Up be the prescribed inclu- 
sion linear displacement, then the force Fp acting on the 
inclusion is 

Fp = —iusZup. (5) 

Fp can be thought of as the resultant of the reactive forces 
from the solid matrix acting on the inclusion, with zero 
net moment. The inverse Z~^ exists since we may con- 
sider the impedance as defined by an imposed force, re- 
sulting in the inclusion displacement Up = (— icjZ)^^Fp. 
We will prove that Z is symmetric (Lemma 1). 

The moment tensor Z relates the moment of the force 
on the inclusion with the inclusion angular velocity, 
—iLuTiBp, 

Mp = -iiuZBp. (6) 
Z is also assumed to be invertible, and will be shown to 
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be symmetric (Lemma 1). 

2. Internal impedances 

The inertial properties of the inclusion are defined by 
two impedance matrices Z p and Z p associated with Un- 
ear and rotational motion, respectively. We call these 
internal impedances since they depend entirely on the 
inclusion and are independent of the exterior region. 

Zp is a mass- like impedance. For a normal solid par- 
ticle it is defined by the mass m as Zp = iajm I. We 
will generally denote Zp as a tensor to include the possi- 
bility of internal structure, although it may be assumed 
on general principles that the impedance is symmetric, 
Zp = Zp. 

Zp is the moment of inertia tensor, and is also symmct- 
~ ~t 

ric Zp = Zp. It has dimensions of a moment of inertia, 
i.e. mass x (length)^. 



Fp= -iwZ Up 




FIG. 3. The inclusion undergoes time harmonic linear dis- 
placement Up, resulting in the net force Fp acting on the 
inclusion. Conversely, if the force Fp is applied then the dis- 
placement is Up. 



3. Summary of main results 

The first principal result is a pair of relations (i) be- 
tween the displacement and force tensors and (ii) between 
the rotation and moment tensors 

Zpi#(r) = (Z + Zp)-iU*(r), (7a) 
Zp^*(s) = (Z + Zp)-iW*(s). (7b) 

Thus, $ = U' and \l/ = W* if the inclusion is immovable 
(infinite impedance), however, the relations (7a) and (7b) 
arc obviously far more general. An immediate corollary 
is that the motion of the inclusion caused by the remote 
force at s is 

Up = {iuj)-^{Z + Zp)-^V\s)F, (8a) 
Op^ {iLj)-\Z + Zp)-^W\s)F. (8b) 



Equations (7a) and (7b) are proved in the next Section 
(Lemma 2). 

The second set of principal results concern applica- 
tions to a spherical inclusion in an isotropic matrix. 
The displacement and rotations tensors U and W have 
particularly simple forms when expressed in terms of 
some lumped parameter impedances introduced in Ref. 6. 
Combined with Eqs. (7a) and (7b) these lead to a se- 
ries of useful identities for the force, moment, displace- 
ment and rotation of the sphere under different excita- 
tion. For instance, the total force and moment on the 
sphere caused by a time harmonic longitudinal or trans- 
verse plane wave is 

Fp = Auo, Mp Fuo A n, (9) 

where Uq is the wave displacement at the center when 
the sphere is not present, n is the propagation direction, 
and the scalars A, F, depend upon the wave frequency, 
particle radius, and other material parameters according 
to Eqs. (44), (46) and (48). Equations (9) could be called 
Faxen relations for solids, by analogy with the use of the 
term in viscous fluid dynamics. 

III. RECIPROCITY BASED IDENTITIES 

Several identities are derived in this Section: (i) the 
symmetry of the external impedance matrices Z and Z, 
(ii) the relation (7a) between the force and displacement 
tensors, $ and U, and (iii) Eq. (7b) relating the moment 
and rotation tensors \I/ and W. The common theme is 
the use of the dynamic reciprocity. 

Consider two distinct fields in V, labeled j = 1, 2, each 
with displacement u^^' , stress crt-') and applied body force 
density per unit volume f '■-'^ all in dynamic equilibrium, 

div cr'^') + ptj^uO) ^ fU) ^ 0, in V. (10) 

The reciprocity identity (Betti's theorem) follow from 
standard arguments^, 

/ dsu(i).r(2)_ / dz;u(i).f(2) 

JdVp Jv 

= f dsu(2).r(i)- / d«u(2).f(i), (^^^ 
JdVp Jv 

where r is the traction vector. The surface integrals in 
Eq. (11) involve only quantities in the matrix. We as- 
sume that the following conditions hold at the interface 
dVp'. (i) continuity of traction, and (ii) the exterior dis- 
placement u is related to the inclusion displacement up 
by 

u - Up = Ar, re dVp, (12) 

where A = A* is a material parameter. This spring- 
like interface condition allows for the possibility of, for 
instance, tangential slip, which we will include in the 
example of the spherical inclusion later. For the moment 
we leave A as arbitrary. 
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The rigid body motion of the inclusion for each of the 
two distinct sohitions in Eq. (11) is assumed to be a 
linear displacement Up-*' and a twist 0^^\ such the to- 
tal displacement is Up^ = Up^ + A r, see Eq. (2). 
Substituting for u'--'-' in the surface integrals then gives 

JaVp JdVp 

- ( di;u(i).f(2) / dz;u(2).f(i) 

Jv Jv 

+ u(2)./' dsT(i)+0f. / dsrArd). (13) 

Note that the interfacial tensor A does not appear in this 
identity. We are now ready to derive the fundamental 
relations, first considering the impedances. 

A. Symmetry of the external impedances Z and Z 

Assume that no force acts in the solid for both solu- 
tions, so that f^^^ = f^^^ = 0. Then Eq. (13) reduces 
to 

u«./ dsT(2)+0W./' dsrAT(2) 

JdVp JdVp 

= 42)./ d5T(i)+0(2). /■ dsrAr^i). (14) 
JaVp JdVp 

The integrals produce the resultant force and moment 
on the inclusion, which follows from the definition of the 
impedances Z and Z as, 

I dsr(j') = -iwZu(/\ (15a) 

JdVp 

I d s r A t(J) = -iujZel^'^ , (15b) 

JdVp 

for j = 1, 2. The reciprocity relation becomes 
uW.Zu(2)+0W.Z0f =u(2).ZuW+0(^).Z0«. (16) 

Since 6^^^^ and 0^^^^ are arbitrary, we deduce 

Lemma 1 The linear and rotational external impedances 
are symmetric, 

Z = Z*, Z = Z*. (17) 

B. Relation between the force and displacement tensors 

We again take field 1 as the solution for the inclu- 
sion undergoing arbitrary rigid body displacement Up^ = 

u^^' -t- e'p^'^ A r with f = 0. Let field 2 be the solution 
for a point force F at s: 

f(2)(x) = F(5(x-s), seV. (18) 



The solution to this, u^^\ is in fact the Green's function 
in the presence of the movable inclusion. Our objective 
is to avoid explicit calculation of the Green's function. 

The displacement on the inclusion surface is again a 
rigid body displacement, Up' ^ Up ' + 6^ ' Ar , and there- 
fore the reciprocity identity (13) becomes 

uW./ dsr(2)+0W./ dsrAT(2)-F.u(i'(s) 

JdVp JdVp 

= u(2)./ dsr«+0(2)./ dsrArW. (19) 

JdVp JdVp 

The integrals involving x^^-* again give resultant force and 
moment according to Eqs. (15) with j ~ I. For field 2, 
let Fp and Mp denote the resultants caused by the point 
force at s, 

[ dsr(2)=Fp, (20a) 

JdVp 

[ dsr Ar(2) = Mp. (20b) 

JdVp 

The displacement at s for field 1 follows from the defi- 
nition of the tensors U and W as u(^)(s) = U(s)up^'' -|- 
W{s)e'j^\ see Eq. (3). 

Elimination of these quantities from Eq. (19) implies 

Fp • uW + Mp ■ 0W - F • U(s)uW - F • W(s)0« 
= -zc^u(,2) . ZuW - iujO'^^ ■ Zel^^ . (21) 

But the rigid body displacement Up^-* and twist are 
arbitrary, and using the symmetry of the impedances, we 
deduce 

Fp = U*(s)F-iwZu(,2), (22a) 

Mp = W* (s)F -iujZe^^l (22b) 

A second set of independent relations follow from the 
equilibrium of the inclusion, or Newton's second law ap- 
plied to a rigid body, 

Fp^iujZpu^^\ (23a) 

Mp^iojZpe^^K (23b) 

Eliminating the linear displacement Up and twist 0^ ' 
between Eqs. (22) and (23) gives 

Fp = Zp(Z + Zp)-iU*(s)F, (24a) 

Mp = Zp(Z + Zp)-iW*(s)F, (24b) 

Finally, referring back to the definition of ^> and Vl/ in 
(4) implies the desired relations: 

Lemma 2 The displacement and force tensors are re- 
lated by 

*(s) = Zp(Z + Zp)-iU*(s). (25) 

The rotation and moment tensors are related by 

*(s) = Zp(Z + Zp)-iW'(s). (26) 

We are now ready to examine these quantities for a 
particular case, the spherical inclusion. 
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IV. SPHERICAL INCLUSION, ISOTROPIC MATRIX 

A. Definition of the problem 

The inclusion has radius a and is embedded in a uni- 
form isotropic elastic medium of infinite extent with mass 
density p and Lame moduli A and p. The interface condi- 
tions at r = a are: (i) continuity of normal displacement, 
(ii) satisfaction of a slip condition. The latter allows for 
relative tangential slip between the inclusion and matrix, 
and is defined by the tangential component of the trac- 
tion T = erf where a is the stress tensor and f = r~^r 
denotes the unit radial vector. The tangential component 
satisfies 



i = zi (vp - v) • t, 



(27) 



where t is a unit tangent vector, v the velocity of the 
elastic medium adjacent to the sphere, vp = —ituup is 
the total velocity of the inclusion at the interface r = a, 
and zi is an interfacial impedance, introduced in Ref. 6. 
This corresponds to A = {iujzi)^^(l — n^n) in Eq. (12), 
where n is the interface normal. The results of Section 
III therefore apply for this slip condition. 

In summary, the conditions at the surface of the sphere 
are 



u • r = Up • r 

T • t = iujzj (u — Up) • t 



(28) 



B. External impedances 

Symmetry arguments imply that the net force (mo- 
ment) exerted on the sphere by the surrounding medium 
and the resulting linear displacement (axis of rota- 
tion) are parallel. Hence, the external impedances are 
isotropic. 



ZI, 



ZI. 



(29) 



The linear impedance Z has been considered previously^ 
while the rotational impedance Z is new. Expressions for 
both are given next. 



1. The linear impedance 

The scalar Z can be expressed in a form reminiscent 
of lumped mass systems^ 



1 



Z + Zm Zl + Zm Zs + Zm 
where the additional impedances are 

7 ■ 4 3 

Zm =iuj-na p, 

Zl = (iw)"M7ra(A + 2p) (1 - ika), 
Zt = {iuj)~^ 4:nap{l — iha), 
1 _ 1 ^ 1 
Zs Zt Atto^zi + (ia;)^^ Snap 



(30) 

(31a) 

(31b) 
(31c) 

(31d) 



Here k and h are, respectively, the longitudinal and 
transverse wavenumbers, k = lu/cl, h = lo/ct with 



cl = \/(A + 2p)/p and ct = ^ pj p. The impedances 
in (31) depend upon and are defined by the matrix prop- 
erties, except for Z5, which involves the interface viscos- 
ity term z/. Thus, Zm is the mass- like impedance of a 
sphere of the matrix material of the same size as the in- 
clusion. Note that Zg = Zt if the inclusion is perfectly 
bonded to the matrix (z/ — > 00). See Ref. 6 for further 
discussion of this and other limits. 



2. Tlie rotational impedance 

The rotational impedance of a spherical inclusion has 
not, to our knowledge, been presented in the literature. 
A derivation is given in Appendix A, with the result that 
Z is 



i(l — iha) 



(32) 



Z Sira^zi Zm + Zt 
The parameters in this identity were defined previously. 

C. Displacement, rotation, force and moment tensors 

1. Internal impedance 

The internal impedances Zp and Zp are necessary in 
order to relate the displacement /rotation tensors with 
the force/moment tensors via Lemma 2. For the sake 
of simplicity we restrict consideration in this paper to 
internal impedances that are isotropic: 

Zp = ZpI, Zp = ZpI. (33) 

For instance, a uniformly solid sphere of mass m has 



Zp = iLum, Zp = iuj-a m. 

5 



(34) 



2. Linear motion 



The displacement tensor U of Eq. (3) is derived in 
Appendix B as 



U(r) = a{Z + Zm) 

^ 1 Zs 

Zs + Zm Zt 



1 



VV- 



Zl + Zm 



VV + h'l) 



hh 



r>a. (35) 



The force tensor $ of Eq. (4) follows from Lemma 2 
and the fact that the impedances satisfy (29) and (33). 
Thus, 



$(r) = y^U(r). 

Z + Zp 



(36) 



The displacement and force tensors satisfy U(— r) = 
U(r), 4'(— r) = *(r) and are symmetric, U = U*, $ = 
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$ . We focus on the properties of U since those of # are 
easily obtained through (36). 
Equation (35) imphes 



where /i„ are spherical Hankel functions of the first kind^ 
and f = r~^r denotes the unit radial vector. In particular 
ho{z) = (iz)~^e*^. In expanded form, 



U = 



Z + Zm r hi{kr) 



Zl + Zm '-krhoika) 
Z + Zm Zs ,,ho{hr) 



/i2(fcr) 
ho{ka) 



f (x) r 



hiQir) 



Zs + Zm Zt ho{ha) hrho{ha) 
^ h2{hr) ~ 
ho{ha) 



(37) 



U(r) = {Z + Z, 



1 



Zl 



Zm 
1 



[( 



Zs + Zm Zt 



1 i 
(kr)^ kr 

^ nr 



)I+(1 



(hr) 



3i 
kr 



(kr) 
(1 



)f ® f] e 



ik{r—a) 



{hr) 



r)f< 



) re 



ih{r—a) 



(38) 



3. Rotational motion 

The skew tensor W relating the rotation to the dis- 
placement at a distance follows from Appendix A as 



hi[na) Hira zi 



(39) 



Then which relates the moment on the inclusion to 
an applied force at a distance, is 



*(r) = ^ W(r). 

Z + Zp 



(40) 



The rotational tensors satisfy are odd functions of their 
arguments, W(— r) = — W(r), '$'(— r) = — 'J'(r), and are 
skew symmetric, W = -W*, * = 



V. APPLICATIONS 

This Section explores implications of the general the- 
ory to the particular case of the spherical inclusion. 



A. Force on a particle from plane wave incidence 

The force on a particle due a remote point load is given 
directly by the tensor <&(r). Taking the source point to 
infinity the effect of the excitation on the particle is equiv- 
alent to an incident plane wave, or a combination of two 
incident plane waves. The far-field form of $(r) follows 
from Eqs. (36) and (38) as 



*(r) =Zp 



Z + Zm \ a 
Z + Zp 

^ih(r-a) 2s ^ 

Zs + Zm Zt 



i 

Zl + Zm 



+ 0(r-2). (41) 



At the same time, the far-field free space Green's function 
is (sec Eq. (61)), 



47r/xr 



-v(g)i) +0(r"^). 

(42) 

Consider, for instance, a unit point force in the far- 
field at r in the direction n = — f . This produces a 
longitudinal plane wave at the origin of the form u = 
u^ne^^^'-^ where uq = e^''^' / {infiK^r). The force on the 
spherical particle due to an incident longitudinal plane 
wave 



u(x) 



ifcn-x „ 
e uo. 



Uo A n = 0, 



is therefore 



(A 4- 2^) 



AnaZp f Z + Zm 



er''"' Uo. 



(43) 



(44) 



Z + Zp \ Zl + Zm 

In the same manner, the force on the spherical particle 
due to an incident transverse plane wave 



u(x) 



e Uo, 



Uo ■ n = 0, 



IS 



AuaZp Zs f Z + Zm \ ^iha 



Z + Zp Zt \ Zs + Zj^i 



(45) 



(46) 



The values of the plane wave induced forces for the rigid 
immovable particle follow from Eqs. (44) and (46) in 
the limit as Zp — * oo. These values actually coincide in 
the static limit, as discussed below after we consider the 
quasistatic limit of Z . 

Nagem and Davis^*^ considered plane wave incidence 
on an elastic sphere in a compressible viscous fluid, with 
specific results focused on the rigid immovable limit. This 
is equivalent to an isotropic clastic medium with shear 
modulus /i = —iLupiyo, where fo is the kinematic viscosity, 
and with a viscously damped longitudinal wave. Their 
expression for the force on the rigid sphere under acoustic 
plane wave incidence (Eqs. (30) and (31) of Ref. 10) 
should agree with (44) in the rigid limit. 
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B. Moment on a particle from a plane wave 

The far-field form of the moment tensor is, from Eqs. 
(39) and (40), 



*(r) 



aih{r-a) ^ 



1 — {iha)~ 



Z + Zt 



axt(f). (47) 



Based on the discussion for the forcing from plane wave 
incidence, it is evident that a longitudinal wave produces 
zero net moment on the spherical particle. A transverse 
plane wave, does however, exert a moment. It may be 
shown that the plane wave (45) produces 



' e ^''''uoAn. (48) 



1 - (i/ia)-i Z + Zp 
The rigid and quasistatic limits are discussed below. 

C. Rigid body displacement due to a plane wave 

The particle displacement under plane wave incidence 
is a combination of a linear displacement and a rigid body 
rotation. These follow, respectively, from the forcing Fp 
of Eqs. (44) or (46) and the moment Mp of Eq. (48) as 

Up = {iiuZpy p + {iujZpy^Mp Ar, r < a. (49) 

Symmetry dictates that the moment tensor Mp is zero 
for longitudinal incidence. 



The expansion (50) goes further than in Ref. 6 (Eq. 
(30)) which did not contain the Cy term. If the low 
frequency expansion is of the form Z = Z^^^"^ {iuj)~^ + 
^(0) 4- Z^^\iuj) + . . . then the coefficient Z'^' determines 
the extra inertia or added mass caused by the linear mo- 
tion of the infinite matrix. The virtual mass coefficient is 
defined as Z'-^'>/Zm, and is therefore Cy of Eq. (53). As 
shown in Fig. 4, the coefficient is positive under no slip 
conditions for all permissible values of Poisson's ratio. It 
approaches the limiting value of Cy ^ 1/2 in the limit of 
incompressibility, j/ — > 1/2, in agreement with the value 
for viscous fluids^^. In contrast, the virtual mass coeffi- 
cient is always negative when the inclusion is permitted 
to slip, and is always less than the incompressible limiting 
value of Cy = —1/6. 




D. Quasistatic limit 

1. Linear motion: Virtual mass 

The quasistatic limit of vanishingly small frequency 
(w — > 0) yields 



127ra/i 



Cy 



1 



{hay + 0(/i^a'') 



iha{2 + K~^) 
(2 + x + «-2)2 

, |/ia|,|A:a| < 1, (50) 



where 



CL 

ct 



2(1-^) 
1 - 2v 



(51) 



and V is the Poisson's ratio. The non-dimensional factor 
X is related to the interface impedance zj in this limit, 
and is chosen so that it takes on the values zero or unity 
in the limit that the sphere is either perfectly bonded or 
perfectly lubricated. 



X 



0, no slip, zi oo, 

1, slip, zj = 0. 



The parameter Cy is 

(2 + K-3)2 



Cy 



i2 + x + K-^y (2-fx + K-2) 



{2^lx + K~^) 

2^2" 



(52) 



(53) 



FIG. 4. The virtual mass coefficient Cv of the spherical inclu- 
sion plotted as a function of the Poisson's ratio. The no-slip 
and slip curves correspond to x = and x = 1 in Eq. (53), 
respectively. 



2. Linear motion: Static displacement 
Based on Eq. (50) we have 



U =- 



3a 



2 + x + i^- 



-I + lim VV-4- I 1 

r h^a h'^r V 1 — iha — ^h^a^ 



1 — in ^ha — '^h'^a^ 



= ^" _j ll+l(.-^-l)VV(r + f)]. (54) 
2 + x + K^r 2 3r 

Evaluating the gradients and removing k in favor of v 
gives 

= 2.[5-6. + 2x(l-^)] ^^' " + ^ 

+ (l 5")^" ® 1 static & immovable. (55) 
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This can be compared with Walpole's^^ result for the 
static perfectly bonded immovable spherical inclusion, 
Eq. (3.21) of Rcf. 12. He considered the force tensor 
which as we have seen is equal to U* in the limit of a 
fixed and immovable rigid inclusion, i.e. Zp ^ oo. But 
U is symmetric for the sphere, and therefore Eq. (55) 
represents both U and When x ~ this agrees with 
Walpole. The static result for x = 1 appears to be new. 

A more precise definition of the fixed inclusion limit is 
that Zp/Z — > oo. At the same time we are taking the 
static limit, so the simultaneous static and immovable 
limit is 



lim 



+0 ujZp 



3. Rotational motion: Virtual mass 



0. 



(56) 



In the limit of low frequency the rotational impedance 
Z of Eq. (32) for approximates as 



Z 



87ra' fj. 

0, 



no slip, 
slip, 



(57) 



where slip and no slip correspond to the limits zj = 
and zj —> oo, respectively. The term iuj^ira^p = 2a^ZM 
can be identified as the virtual mass due to the rotating 
solid. The internal rotational impedance of a solid sphere 
is Zp = ^o? Zp. Hence the virtual mass in rotation is five 
times the mass of solid matrix material in the volume of 
the sphere, in agreement with a similar result for Stokes 
flow". 



4. Quasistatic plane wave force on an immovable particle 

The forcing on the particle from plane wave incidence is 
the same, whether longitudinal or transverse waves are 
incident, in the quasistatic limit for a rigid immovable 
sphere. It may be checked that both (44) and (46) be- 
come 



127ra/i 



X + K 



Uq. 



(58) 



This apparently strange result may be reconciled with 
the physical nature of the limit: the matrix moves by the 
static displacement Uq, while the particle is stationary. 
It is therefore sensible that there is a net force on the 
particle, and that it is in the direction of Uq. This also 
explains the identical form of the limit for both types of 
wave incidence. The static limit is one of the unusual 
features of the the immovable sphere. For further discus- 
sion sec Pao^'^, who quite properly questions the physical 
validity of the rigid fixed assumption. Among its failings, 
as Pao notes, this configuration does not display Rayleigh 
scattering behavior at low frequencies. 



5. Quasistatic plane wave moment on a fixed sphere 

The moment on the rigid immovable spherical particle 
has a particularly simple form, 

= — !^47ra2u n e"'^"- uq A n, immovable. (59) 
1 — iha 

This is valid at all frequencies, but as a; ^ it vanishes, 
unlike the force on the particle in the same limit. 



E. Small inclusion limit 

1. Linear motion 

As a ^ we have Z^, Zg, Zt = 0(a), while Zj^j, 
Zp = o{a). It may be easily verified that U and $ reduce 
in this limit to 

U(r) = iwZG(")(r), ^{r) ^ iujZpG^^\v), (60) 
where G'"-' is the free space Green's tensor 



G(")(r 



1 

A'Kp 



^ihr 



-1 + vv 



Akr ^ 



(61) 



In hindsight, the form of $ is obvious based on the dy- 
namic equilibrium of the inclusion : iuZpVLp 
Up = G("'F, Fp = *F, and U then follows from (36) 



Fp with 



2. Rotational motion 



The rotational quantities W and on the other hand, 
become negligible in the small inclusion limit. This fol- 
lows from the scaling W =0(a) in Eq. (39). 



F. Surface displacement and traction 

1. Linear motion 

The displacement tensor U, which is defined in the 
exterior region, reduces to the following on the interface: 



U(af) = 1+ ( 



Zs — Zt ^, Z + Z]s,i 
Zt Zs + Zm 



)(l-f®f). (62) 



This becomes the identity I under no-slip conditions, 
since then Zg — Zt ~ 0. Alternatively, the interface 
conditions (28) can be written 

u = Up + (icc)Z/)^"^(-t (g) t)T, r = a. (63) 

The final term on the RHS vanishes as z/ oo, which is 
the no-slip limit. 

Substituting u = U(af)up on r = a, in (63) provides 
an explicit expression for the interfacial shear traction in 
terms of the linear displacement Up , 

,Zs — Zt-^^ Z + Zm 



T ■ t 



VjJZi 



(64) 



Zt Zs + Zm 

The shear traction vanishes under pure-slip conditions 
{zi = 0), and for a bonded interface it becomes 

Zt / + Zjyi , 



(65) 



Faxen relations in solids 8 



2. Rotational motion 



2. Spherical inclusion 



The displacement on r = a is, 

W(af) = —afl axt(f), 



(66) 



where fJ, given in (39), reduces to unity under no-slip 
conditions. Conversely, f2 = for pure slip, indicating 
that the solid does not move even as the inclusion rotates. 
In this case the traction is pure shear, and 



. 3Z ^ . 

-^tijJ- — T 0„ A r, 



(67) 



This is non-zero except under pure slip conditions, when 
Z ^ 0, and there is no rotational interaction between 
the inclusion and the matrix. 



G. Acoustic limit 

1. General formulation 

Finally, we discuss how the general elastodynamic for- 
mulation reduces when the matrix is an acoustic fluid. 
In this limit shear effects are ignorable and the medium 
is characterized by density p and bulk modulus K = pc^, 
where c is the acoustic wave speed. Taking the displace- 
ment u and pressure p as field variables, the momentum 
balance and constitutive law are, respectively, 



uj pu~ Vp, p = — A'V • u. 



(68) 



The acoustic wavenumber is k = loj c. 

We introduce two vector functions q(r) and (/)(s) that 
are analogous to the tensors U and If the inclusion 
is moved back and forth with the displacement Up then 
the condition on the inclusion surface is that the normal 
velocity is continuous. 



u • n = Up • n on dVp. 



(69) 



The pressure at a point r away from inclusion is defined 
by q as 



p{v) = q(r) • Up rin V. 
Conversely, consider a voluminal source at s: 

V2p + fc^p = /(5(x-s). 
The force on the inclusion is 



Fp = - y" dspn = /</>(s). 



(70) 



(71) 



(72) 



which defines the vector function ^. 

The connection between q and is given by 

Lemma 3 The acoustic displacement and force vectors 
are related by 



0(s) = (p^2)-i2^(2 + Zp)-iq(s). 



(73) 



This may be derived by application of reciprocity to the 
acoustic (Helmholtz) equation, in a manner similar to 
how we derived Lemma 2. 



Finally, we consider the example of the spherical inclu- 
sion. The vector q follows from the acoustic limit of the 
elastic result in (35), 



q(r) = aK\ 



Z + Zm ^ ^ e 



ik{r — a) 



-, r>a, (74) 



where Za, analogous to the longitudinal impedance 
in elasticity, is 



Za = {i^jj) ^4:TraK (1 — ifca). 



(75) 



Zm is as before, and the sphere impedance Z is now given 
by (30) with Zs = 0, which implies 



1 _ 2 ^ 3 
Z Zm Za 



(76) 



VI. CONCLUSION 



Starting from the notion of an inclusion with the six 
degrees of freedom of a rigid body, we introduced dis- 
placement/rotation and force/moment tensors relating 
the motion of the inclusion to the displacement and force 
at arbitrary exterior points. These can be considered 
as generalized Green's functions appropriate to the con- 
strained nature of the inclusion. The general relations 
(7a) and (7b) between the displacement/rotation and 
force/moment tensors are one of the main contributions 
of the paper. These identities are extremely useful in pro- 
viding a means by which one can consider the dynamic 
properties of particles embedded in a solid matrix. 

Useful results have been obtained for the simplest but 
important configuration of a uniform spherical particle. 
The linear and rotational impedances. Z and Z, arc given 
in Sections IV.B.l and IV. B. 2, respectively, the latter for 
the first time. Explicit expressions are given in Eqs. (36) 
- (40) for the displacement tensors U and W and for the 
force tensors $ and ^. Perhaps the most practical new 
results are Eqs. (44), (46) and (48) which provide simple 
formulae for the force and moment on a particle under 
plane wave incidence. The associated displacement of 
the particle is given by Eq. (49). These concise expres- 
sions resemble Faxen relations that are frequently used 
in microhydrodynamics. 

Equation (50) extends the quasistatic expansion of 
Norris^ to include the virtual mass coefficient, which can 
be negative if slip occurs. The quasistatic form of U, Eq. 
(55), which relates the displacement of the inclusion to 
particle displacement in the matrix, generalizes a recent 
formula of Walpole^^. The quasistatic limit of the plane 
wave force on a sphere, Eq. (58), is reminiscent of Stokes 
drag law, but is proportional to the displacement vector 
of the incident plane wave. It also includes the possibility 
of slip relative to the matrix. However, the moment on 
a rigid sphere from plane wave incidence is proportional 
to the incident particle velocity, Eq. (59), and vanishes 
in the limit of zero frequency. Other limiting cases con- 
sidered include the small inclusion limit, and the purely 
acoustic limit. 
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Taken together the resuhs of paper offer a consistent 
means for analyzing wave-particle interaction in elastic- 
ity and viscoelasticity. Future applications will look at 
replacing the solid spherical particle with more compli- 
cated, and more interesting, internal structure. This 
amounts to considering more general forms of the inter- 
nal impedances. The results developed here can also be 
used to develop simplified methods for scattering from 
particles. These issues will be examined in separate pa- 
pers. 



APPENDIX A: ROTATION OF A SPHERE 

The sphere r < a undergoes oscillatory rotation up = 
9p A r. Let e be the axis of rotation, so that 9p = 6pe, 
and consider the possible solution u = V A e/ in the 
matrix r > a. This satisfies the equations of motion 



u + VV ■ u - /i"^ V A V A u = 0, 



(Al) 



provided that / is a solution of the reduced Helmholtz 
equation V^/ -t- h^f = 0. The function / must depend 
only on r in order to match the prescribed rotation on 
r ~ a. Hence, 

u = phi Qir) e A f , r > a, (A2) 
where r = r/r. The traction in an isotropic solid is*", 

d 1 

T = fAdivuH — gradr-u + /i(— )u, (A3) 

r or r 

from which we obtain 

3 

T = l3^h\ho{ha) ~ —hi(ha)\e ^Y. (A4) 
ha 

The boundary conditions (28) therefore reduce to a single 
equation for the parameter (3 of Eq. (A2), 



/3^h\ho(ha) — - — hi(ha)] = iuzjlphilha) ~ aOp\. (A5) 
ha 

The moment Mp = Jdsr A t is obtained from the 
identity 



f 8 

/ dsr A (e A f) = -TTo^e, 



(A6) 



and the impedance then follows from the definition (6) 
as Z = where Z is given by (32). 

APPENDIX B: DISPLACEMENT OF A SPHERE 

The solution to the radiation boundary value problem 
of the spherical inclusion undergoing linear displacement 
with boundary conditions defined by (28) has been solved 
by Oestreicher^ for the case of no slip and more recently, 
by Norris^, with slip included. We follow the latter with 
some slight changes in notation. The solution is based on 



the following representation for the elastic field outside 
the sphere, r > a, 

u = -Ci Up • yVhoikr) + C2 Up • (VV - IV'^)ho{hr), 

(Bl) 

where r ~ |r| is the spherical radius and ho is the spher- 
ical Hankel function of the first kind, ho{z) ~ (iz)~^e'^. 

In the notation of Ref. 6, Ci = — Ai/fc^ and C2 = 
3Bi/h^. Also, with reference to Norris^, the inclusion 
displacement is of magnitude uq: Up = mqX. Equations 
(15) and (19a) of Ref. 6 combined with Eq. (29), give 
(noting that of Ref. 6 is now Zm ) 



h2ika)Ai = {l + —)uo, (B2a) 

6^Bi = - -^uo. (B2b) 

ha ka Zm 

Using hi{z) = — /iq(z) and /i2 = — /iq + 3z^^hi implies 
the identities 



hi{ka) 

kaho{ka) 3Zm' h(){ka) 



Zl h2{ka) Zl + Z^ 



M 



(B3) 



Similar identities for arguments ha instead of ka have Zt 
instead of Z^. Combining these results implies 



Ci = 
C2 = 



z + z_ 



M 



1 



Zl + Zm k'^ho{ka) ' 
Z + Zm n Zs 1 



Zs + Zm Zt h'^ho{ha) 
The expression (35) for U then follows. 



(B4a) 
(B4b) 
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